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Abstract 

In this paper we make a comprehensive study about the four B — t- Kr]^'^ decays in the pertur- 
bative QCD (pQCD) factorization approach. We calculate the CP-averaged branching ratios and 
CP-violating asymmetries of 5 — t- Kr]^'^ decays in the ordinary rj-r]' mixing scheme, the rj-rj'-G 
and the r]-rj'-G-gc mixing scheme, respectively. Besides the full leading order (LO) contributions, 
all currently known next-to-leading order (NLO) contributions to i? — )• Kr]^'^ decays in the pQCD 
approach are taken into account. From our calculations and phenomenological analysis, we find 
that (a) the NLO contributions in general can provide significant enhancements to the LO pQCD 
predictions for the decay rates of the two B — )• Kr]' decays, around 70% — 89% in magnitude, but 
result in relatively small changes to Br{B — )• Kg); (b) although the NLO pQCD predictions in 
all three considered mixing schemes agree well with the data within one standard deviation, but 
those pQCD predictions in the g-g'-G mixing scheme provide the best interpretation for the mea- 
sured pattern of Br{B Kg^'^): Br{B^ K^g) ^ 1.13 x 10-^ Br{B^ K^g') ^ 66.5 x 10"^, 
Br{B^ K^g) ^ 2.36 x 10"^ and Br{B^ K^g') » 67.3 x 10"*^, which agree perfectly 
with the measured values; (c) the NLO pQCD predictions for the CP-violating asymmetries 
for the four considered decays are also in good consistent with the data; (d) the newly known 
NLO contribution to the relevant form factors Aipp can produce about 20% enhancement to 
the branching ratios Br{B — )• Kg'), which plays an important role in closing the gap between 
the pQCD predictions and the relevant data; (e) the general expectations about the relative 
strength of the LO and NLO contributions from different sources are examined and confirmed 
by explicit numerical calculations. 
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I. INTRODUCTION 



From the first observation of unexpectedly large branching ratios for B — Kr]' decays 
was reported by CLEO in 1997 Q, the four B — t- Kr]^'^ decays have been a "hot" topic for 
a long time. Although many physicists have made great efforts to explain the pattern of 
very large branching ratios Br{B — )■ Kr]') and very small branching ratios Br{B — )■ Krj) 
[2|-l9|, but it is still a hard job to understand such pattern in the framework of the standard 
model(SM). 

On the experimental side, the branching ratios of the four B — )■ Ki]^'^ decays have been 
measured with high precision (lol . 11|: 



5r(5° ^ /^%) = (1.23l°;2^) x 10-^ 
fir(5° ^ K%') = (66.1 ± 3.1) X 10"^ 
BriB^ ^ K^T]) = i2.3Qtoil) x 10-^ 

Br{B^ K^T]') = (71.1 ± 2.6) x 10-^ (1) 

For the relevant CP-violating asymmetries, currently known experimental measurements 
are 0, [ll| 

A'^'piB'' -> irV) = (1 ± 9)%, 
^™^(5° -> i^V) = (64 ± 11)%, 
A'^jpiB^ ^ K^V) = (-37 ± 8)%, 

A'SUB^ ^ K^V') = (l-3l}:?)%, (2) 



On the theory side, these decays are calculated recently in Ref. [121] by employing the 
perturbative QCD factorization approachf9'] and using the ordinary Feldmann-Kroll-Stech 
(FKS) rj — r]' mixing scheme with the inclusion of the partial NLO contributions, i.e., 
the QCp vertex corrections, the quark- loops and the chromo- magnetic penguins Ogg- In 
Ref. [12| the authors found that the NLO contributions can provide a 70% enhancement 
to the LO pQCD predictions for Br{B — )■ Kr]'), but a 30% reduction to the LO pQCD 
predictions for Br{B — )■ Kr]) [ll]: numerically Br{B^ — )■ K^r]) f» 2.1 x 10~^, Br{B^ — )■ 
K'^r]') ^ 50.3 X 10-^ Br{B+ K+r]) ^ 3.2 x 10"^ and Br{B^ K^r]') ^ 51.0 x 10^^; 
although the differences between the pQCD predictions and the data decreased effectively 
by the inclusion of the partial NLO contributions, but the central values of the pQCD 
predictions for Br{B — )■ Kr]') are still lower than the data by about 30%. As for the 
CP-violating asymmetries, the pQCD predictions in Ref. [l2[ are basically agree well with 
the data. 

Very recently, three new progresses in the studies of the two-body charmless hadronic 
B — >■ M2M3 decays (here Mj stands for the light mesons, such as n, K,r]^'\ etc.) in the 
pQCD approach have been made: 



In Ref. [ij], Li et al. calculated the NLO contributions to the form factors of S — vr 
transitions in the pQCD approach and found that the NLO part can provide about 
20% enhancement to the LO one. The enhanced form factors may then lead to 
a larger branching ratio for B — )■ Kr]' decays. The still missing NLO parts in 
the pQCD approach are the O(a^) contributions from non-factorizable spectator 
diagrams and the annihilation diagrams, they are most possibly small according to 
general arguments. 
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(ii) In Ref. [15|, the authors studied and provided a successful pQCD interpreta- 
tion for the Belle measurements of Bd/Bg — > J/'i!T]^'\ i.e., R^^'^ = Br{Bq — > 
J I'^r]') / BriBq — )• J/'^rf) < 1 with q = {d,s), by using the f] — rj' — G mixing 
formalism [l6|], where G represents the pseudoscalar glueball. This result encourage 
us to check the possible effects of the pseudoscalar G on _B — )• Ki]^'^ decays, although 
such contributions may be small as generally expected |17l] . 



(iii) In Ref. [18|], the authors studied the rj — rj' — G — rjc mixing scheme, obtained 
constraints on the mixing angle (pQ, (pc ~ 11°, between G and rjc by fitting to the 
observed rjc decay widths and other relevant data, and found that the t]c mixing can 
enhance the pQCD predictions for Br{B — )■ Krj') by 18%, but does not alter those 
for Br{B Kt]). In Ref. [l8|, the authors superposed the contributions from B — )■ 



K'r]c due to the rjc mixing onto the partial pQCD predictions as given in Ref. [12 
directly. They did not consider the effects of newly known NLO contributions to 
the corresponding form factors F(f^^(0) and F(f~^^*'(0). 

Motivated by above new progresses [3, 1^, 18|, we think that it is time for us to make a 
comprehensive study for the four B — )• Krf^'^ decays in the pQCD factorization approach. 
We will focus on the following points: 



(i) Besides those NLO contributions already considered in Ref. [12[, we here will firstly 
extend the calculation of the NLO part of the form factors for i? — )■ vr transition 
to the cases for the similar B ^ K and B — t- (r/^, rjs) transitions, and then take 
these newly known form factors at NLO level into the calculations for the branching 
ratios and CP-violating asymmetries of 5 — )■ Krj^'^ decays, to check its effects on 
the corresponding pQCD predictions. 



(ii) Besides the ordinary rj — rj' mixing scheme [13|, we will also calculate these four 
decays in the 7]—r]' — G [l6| mixing scheme and the rj—ri' — G—rjc mixing scheme 18|. 
respectively. We want to check the effects of the possible "pseudoscalar Gluball" 
and rjc component of rf^'^ meson on the pQCD predictions. 

(iii) We will calculate numerically for the individual decay amplitudes Ai""^^ ( the emis- 
sion diagrams Figs.l(a)-l(b)), M.^^'^ (the spectator diagrams Figs.l(c)-l(d)) and 
j^anm (^annihilation diagrams Figs. 1(e)- 1(h) ) and compare the relative strength of 
the contributions from different set of the Feynman diagrams at the leading order, 
or from the different sources at the next-to-leading order: such as M.vc (i-e. NLO 
vertex corrections Figs.3(a)-3(d) ) and A^g/ ( i.e. NLO chromo-magnetic penguins 
Figs.3(g)-3(h)). We try to find the source of the dominant contribution, to evalu- 
ate the possible strength of the two still missing NLO contributions in the pQCD 
approach. 

The paper is organized as follows. In Sec. II, we give a brief discussion for the pQCD 
factorization approach and the three different kinds of mixing schemes: rj-rj' , rj-ri'-G, 
and rj-ri'-G-Tjc mixing schemes. In Sec. Ill, we make the analytic calculations of the 
relevant Feynman diagrams and present the various decay amplitudes for the studied 
decay modes in leading-order. In Sec. IV, all currently known NLO contributions in 
the pQCD approach will be investigated. In Sec. V, we will show the numerical results 
for the pQCD predictions for the branching ratios and CP-violating asymmetries of the 
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considered decay modes, we here also calculate and compare the relative strength of the 
LO and NLO contributions from different set of the Feynman diagrams or from different 
sources. The summary and some discussions are included in the final section. 



II. 



THEORETICAL FRAMEWORK 



As is well-known, the pQCD factorization approach has been widely used in studies for 
the two body charmless hadronic B/Bg/Bc — M2M3 decays (here Mi stands for the light 
pseudo-scalar meson P, the vector meson V, the scalar meson S, etc.) js], 0, 0, 12, 14, 19- 
25| . Some pQCD predictions, for example the large CP- violating asymmetries Aqp{B^ — )• 
(30 ± 10)% in Ref. jij and large branching ratio Br{Bs — )■ vr+vr") 5 x 10"^ 



TT ' TT 



Refs. [2 
the four B 
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m 



have been confirmed by experiments [25|]. We here focus on the study for 
Kr]^'^ decays. For more details of the formalism of the pQCD factorization 



approach itself, one can see for example Refs. [ol, [H, 21, 27 



A. Outline of — K-q'^'^ decays in the pQCD approach 

In the B-rest frame, we assume that the light final state meson M2 and M3 is moving 
along the direction oi n = (1, 0, Ot) and v = (0, 1, Or), respectively. We use Xi to denote 
the momentum fraction of the anti-quark in each meson, fcj^ the corresponding transverse 
momentum. Using the light-cone coordinates the B meson momentum and the two 
final state meson's momenta P2 and P3 (for M2 and M3 respectively) can be written as 

i^B = ^(l,l,OT), P2 = ^{l-rlrlOT), Ps = ^{rll-rlO^), (3) 
where rj = rrii/MB with rrii is the mass of meson Mj. If we choose 

/si = ^ [xi, 0, kix) , ^2 = ^ (3^2(1 - rl),X2rl, ksx) , 

h = ^{x,rlxs{l-rl),k,T). (4) 

The integration over the small components fcf , fc^, and will lead to the decay ampli- 
tudes conceptually 



A{B — > M2M3) ~ J dxidx2dxsbidbib2db2b^db^ 

■Tr [C{t)<l>B{xubi)<l>M,{x2,b2)<l>M:Ax3,b3)H{x,,b,,t)St{x,)e-'''^'^] , (5) 

where bi is the conjugate space coordinate of kix- In above equation, C(t) is the Wilson 
coefficient evaluated at scale t, ^b{xi, bi) is the wave function of the B meson, ^m2{x2, ^2) 
and $4/3(0:3, 63) are the wave functions of the final state meson M2 and M3, the function 
St{xi) and e"'^*^*-* are the Sudakov form factors which suppresses the soft dynamics effec- 
tively ja]. 
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For the studied B — )■ Kr]^'^ decays, the corresponding weak effective Hamiltonian can 
be written as 28 : 



K 
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VubV*g 



CMO^M + C2ifi)o^M - Vt,v:AJ2c,if^mf^) +H.c.,(6) 



where q = d,s, Gp = 1.16639 x 10 ^^GeV ^ is the Fermi constant, Oi {i = 1, 10) are the 
local four-quark operators 28|] . The Wilson coefficients Cj in EqEland the corresponding 



renormalization group evolution matrix are known currently at LO and NLO level [28!]. 

In this paper, we will calculate B — )■ Kt]^'^ decays in the pQCD approach with the 
inclusion of all known NLO contributions, focus on the effects of the newly known NLO 



contributions to the form factors Fg^^^(O) and Fq~^^ 



(') 



(0)11 



B. Different Mixing Schemes 

In order to check the effects of different mixing schemes we will calculate the B — )■ Ki]^'^ 
decays in the following three different mixing schemes(MS): 



(i) MS-1: The conventional FKS scheme [13| oirj — rj' mixing in the quark-flavor basis: 
r]q = [uu -\- dd) I \/2 and rfs = ss; 



(ii) MS-2: The rj-rj'-G mixing scheme defined in Ref. [16 

(iii) MS-3: The rj-ri'-G-Tjc mixing scheme defined in Ref. [ll 
Firstly, in the conventional rj — 1]' mixing scheme, the physical r] and rj' can be written 



as 

(7) 



^ 1 = ( + + ^2(0) ss 

V' J ~ [ Fii(l)) iuu + dd) + F^{(f)) ss 



where is the mixing angle, and 

%/2Fi(0) = F^(0)=cos(0), 
F2(0) = -v^F{(0) = -sin(0), (8) 

are the mixing parameters. The relation between the decay constants (/^, /^,, f^,) and 
{fqi fsi) can be found in Ref. 12|. The chiral enhancement itlq and ttIq have been defined 
in Ref. [2l[ by assuming the exact isospin symmetry = = m^. The three input 
parameters /g, fs and in the FKS mixing scheme have been extracted from the data of 



the relevant exclusive processes pi 

/g = (1.07 ±0.02)/,, /, = (1.34 ±0.06)/,, = 39.3° ± 1.0°, (9) 
with /, = 0.13 GeV. 

In the Tj-rj'-G mixing scheme [3], however, the physical states {ri,ri',G) are related to 
{rig,ris,g) through the rotation U{(I),(I)g) 



f/(0,0G) 1^.) , (10) 





with the rotation matrix 



cos + sin ^ sin Ag — sin<^ + sin0 cos^jAg — sin 9 sin (f)G 
sin — cos sin 0j Afj cos0 — cos6' cos6'jAG cos 6* sin 0^ 
— sin 9i sin (pQ — cos 6i sin (pc cos (pc 



(11) 



where AG = 1 — cos0G) = 54.7° is the ideal mixing angle, (t> = 6 + 6i with the 6 varying 
in the range of — 17° < < — 11°. Based on recent analysis |16l. [l8|, we will choose 

^ = -11°, 0G = 12°, fg = U, fs = l.SU, mG = 1.376GeV, (12) 

in the numerical calculations. 

The decay constants associated with the physical {rj, i]', G) states are related to those 
associated with {riq,r]s,g) states via the same mixing matrix 





(13) 



Following Ref. [16|, the chiral enhancement ml and ttiq in the rj-rj'-G mixing scheme 
can be written as 



ml 



mn 



qq 



2m„ 



2m, 
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2m, 2m, 



V2fs 

fq 



Uu), 



22 



V2l 



(14) 
(15) 



with 



Uii = (cos + sin 6* sin OiAcYm'^ + (sin — cos 9 sin 9iAG)'^m'^, + (sin 9i sin 0g')^'^g(?16) 
f^i2 = f^2i = (cos + sin 6* sin 9iA.G) ■ (— sin + sin cos 9iA.Q)m'^ 

+ (sin (j) — cos 9 sin 9iAG) ■ (cos (p — cos 9 cos 9iAG)m^, 

+ sin 9i sin <pG ■ cos 9i sin (pcrriQ , 
f^22 = (— sin + sin 6* cos 9iAGym'^ + (cos (p — cos 9 cos 9iAGYm^, 

+ {cos 9 i sin (pG^mQ . 



(17) 



In the third rj-ri'-G-rjc mixing scheme, finally, the fiavor states are transformed into the 
physical states through the mixing matrix U{9, (pG, > 



( \V) \ 

\G) 

V \Vc) J 



U{9,cPgAq) 



( \Vq) \ 
\Vs) 

\9) 
\ \Vq) I 



(19) 



with the 4x4 mixing matrix 
U{9,<Pg^ 



( c9c9i — s9c(pGs9i 
s9c9i + c9c(pGs9i 

-S(pGs9i 





-c9s9i — s9c<pGc9i —s9s(, 

-s9s9i + c9c(pGc9i c9s(pGC(pQ 

-s(pGc9i cq 

-sc 



-s9s(pGS(pQ \ 
c9s(pGS(pQ 

o 

o 
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where c0 (s0) is the shorthand notation for cos0 (sin0) and similarly for others. The 
decay constants associated with the rj, r]', G, and rjc physical states are related to those 
associated with the rjg, rjg, g, rjQ states through the mixing matrix U{9, (j)G, 

^ fri frj fri\ ( fq fq fq ^ 



fQ fs £c 

fQ fs fc 

Jg Jg Jg 

\ f" f I 



fq f fc 

Js Js Js 

fq fs fc 

Jg Jg Jg 

\f'c fc fcj 



(21) 



Furthermore, we find that the chiral enhancement ml and in the third mixing scheme 
are identical to those as given in Eqs. fll4lll5p . In rj-ri'-G-rjc mixing scheme, we take 
(j)Q = 11° and choose the same input parameters {0, (pc, fq, fs, ^^g) as given in Eq.( IT2|) . 



Wave Functions 



The B meson is treated as a very good heavy-light system. Its wave function can be 
written as the form of 

i 



B 



V2K 



^B + "^B)750B(kl). 



(22) 



Here we adopted the widely used B-meson distribution amplitude in the pQCD approach 



')b{x, h) = Nbx'^{1 — x)^exp 



(23) 



where the normalization factor Nb depends on the value of Wf, and and defined through 
the normalization relation 



dx 05 (x, 6 = 0) 



B 



For the shape parameter Uf, we take Uf, = 0.4 ± 0.04 GeV in numerical calculations. 
The wave functions of the final state mesons M = {K, rjq, rjs) are defined as: 



(24) 



=75 [^20M,(^i) + moi(f)MS^i) - mo^(# - l)(j)lj^{xi)] , (25) 



where rriQi is the corresponding meson chiral mass. Pi and Xi are the momentum and 
the momentum fraction of Mi, respectively. The explicit expressions of the distribution 



amplitudes 



lA,P,Ti 



Xi) for M = {K,riq,ris) will be given at Appendix A. 



In the third r]-r]'-G-r]c mixing scheme, the r/c part of the i]^'^ meson will contribute to 
the B Kr]^'^ decays through the decay chain B — )■ Kr]c — ?■ Kt]^'^ . The wave function of 
the Tjc can be written asjloj: 



=75[^2C(a;2) + m^>^Jx2)]. 



(26) 



c 



The distribution amplitudes (j)'^'^ are of the form |29|: 



zx(l — X 



x(l — x) 



<PLix) = 1.97 



2V2iV; 

fvc 

2v/2iVcLl-2.8x(l-x)J 



0.7 



Ll -2.8x(l -x)J 
x(l-x) i0-7 



(27) 
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(a) (b) (c) (d) 

(e) (f) (g) (h) 



FIG. 1. Feynman diagrams which may contribute to the B — )■ Kr]^'^ decays in the pQCD 
approach at leading order. 

III. B AND B K7]c DECAYS AT LEADING ORDER 

In this section we will present the total decay amplitudes for B — t- Kt]^'^ and B — )■ Krjc 
decays in the pQCD approach at leading order. 



A. B ^ Krj^'^ decays at leading order 

The B — )■ Kt]^'^ decays have been studied by using the ordinary "q-r]' mixing scheme 
and by emp loying the pQCD factorization approach at the LO and partial NLO level 
in Ref. Il2|. We recalculated and confirmed the relevant analytical formulas as given 
in Ref. [l2|. For the sake of the reader, we here present directly the decay amplitudes 
obtained by evaluating the Feynman diagrams Fig.l(a)-l(h). 

For the factorizable emission diagrams Figs.l(a)-l(b), the decay amplitudes for the 
cases of B ^ K transition are of the form: 



:pV-A 
eK 



_pV+A 



dxidxs / bidbibsdbs(j)B{xi,bi) 



X I [(1 + X3)0^(X3) + r3(l - 2X3)(0^(X3) + ^^(xs)) 

xEe{ta)KiXi, X3, 61, 63) + 2r3(f)^{x3)Ee{Qhe{x3, Xi, 63, 61 
/•l /»oo 

F^^ = —lQTir2CFM% I dxidxs / bidbib3db3(f)B{xi,bi) 

Jo Jo 

X { 4>i{x3) + r3(2 + X3)0^(X3) - r3X30^(x3) 
xEe{ta)he{Xi, X3, 61, 63) + 2r30^(a;3)Ee(4)/ie(x3, Xi, 63, 



f28) 



(29) 



where Cp = 4/3, r2 



rriQ^ /Mb and = rx = tuq /Mb, where ttiq* and ttiq are the 



chiral scale parameters. The functions Ee{t), the definitions of the factorization scales t 
and the hard functions are given in Appendix B. 



8 



For the non-factorizable emission diagrams Figs.l(c)-l(d), whose contributions are 

—32'kCfM^/\/6 / dxidx2dxz / bidbib2db2(pB{xi,hi)4>'^{x2) 
Jo Jo 

X I (1 - X2)0i(x3) - r3X3{(f)^{x3) - 4>Kix3)) Kih)h„{xi, X2, X3, hi, 62) 

+ - {x2 + xs)(pi{xs) + r3X3(0^(a;3) + 0^(0:3)) 

xE'Mhn{xi, X2, X3, bi, 62)}, (30) 

rl poo 

— 3271 CpM^r 2/ \/Q / dxidx2dx3 / bidbib2db2(j)B{xi,bi) 
Jo Jo 

X{[(l-X2)0i(x3) «(X2)+<(X2)) 

+^32^3 {(t>^ix2) - (1>J,{X2)) {(t>Ki^3) + 0^(2^3)) 

+ (1 - a;2)r3 (0f (2:2) + 02 (^^2)) (^f (a^s) - E[Xh)K{xi, X2, x^, 61, 62) 

- X20i(x3)«(x2) - <(X2)) + r3X2«(x2) - ^^,{X2)) 

{4>kM - 0^(2:3)) + r3a;3«(a;2) + < (0:2)) (0^(0:3) + (fAx^)) 
XE'^{Qhn{xi, X2, X3, bi, 62) |, (31) 

/•l /"OO 

— 327rCi?M^/\/6 / dxidx2dxz I bidbib2db24>B{xi,bi)4>^{x2) 
Jo Jo 

x| (X2 - X3 - 1)0^(X3) + r3X3(0^(x3) + 0^(X3)) E'^{tb)hn{xi,X2,X3,bi,b2) 

+ X2(f)i{x3) + r3X3{(t)J.{x3) - (t)^{x3)) Eg(4)/i„(a;i, 0:2, 0:3, 61, 62) } , (32) 

where Xi = 1 — Xi, while denotes or 0^^. 

For the factorizable annihilation diagrams Fig.l(e)-l(f), the decay amplitudes are of 
the form 

pi poo 
= FaK"" = -87rCf / dX2dX3 / b2db2b3db3 



M: 



SP 
eK 



X 



(^3 - l)0;^(a;2)0^(x3) -4r2r30j'(a;2)0x(a:3) 

+2r2r3X30j'(x2)(0^(x3) - 0^(2:^3)) Ea{tc)ha{x2,X3,b2,b3) 



+ 



(33) 



a;20;^(x2) 0:^(^3) + 2r2r3(0j'(x2) - 0j(a:2))0K(a:3) 
+2r2r3a;2(0j'(a;2) + 0j(a;2))0^(a;3) -Ea(^c)/ia(^3, 2:2, &3, &2)}- 

/■I /-oo 

Ff/ = -IOttCfM^ / 0^x2(^x3 / b2db2b3db3 
Ja Ja 

2r2(p^ {x2)(f>K{x3) + (1 -a;3)r30;^(x2)(0^(a;3) + 0^(2:3)) Ea{tc)ha{x2,X3,b2M) 
2r30;^(x2)0^(x3) +r2a;2(0j'(a;2) -0j(a;2))0^(a;3)]Fa(tc)^a(a;3,a;2,&3,&2)}. (34) 



X 

+ 



{ 



For the non-factorizable annihilation diagrams Figs.l(g)-l(h), whose contributions are 



M. 



V-A 
aK 



<pBsixi,bi) 



pi poo 

—?)2'kCfM^/\/Q I dxidx2dx3 / bidb2b2db- 
Jo Jo 

X I - X24>^{x2)(f)i{x3) - 4r2r30^(x2)<^^(a;3) 
+r2r3(l - X2)«(a;2) + 0j(x2))(0^(x3) - ^^(xs)) 

+r-2r3a;3(0^(x2) - 0^(x2))(0x(a;3) +0^(3^3)) E'^(tc)hna{Xi,X2,X3,bi,h 
+ X30;J(X2)0^(X3) + X3r2r3(0^(x2) + ^^(Xs)) (^^(Xg) - ^^(Xg)) 
+X2r2r3(0^J'(x2) - 0j(x2))(0x(a;3) + (PxiXs)) 
xK{t'c)h'^a{Xl,X2,X3,bi,b2] 

—32ttCfM^JVQ / dxidx2dx3 I bidbib2db2(f)B{xi,bi) 



(35) 



r2(2 - X2)(0^'(X2) + 0eta^(x2))0^(x3) 



-r3(l + X3)0;J(X2)(0]^(X3) -0^(X3)) E'^{td)hna{Xl, X2, X3, bi, b 



+ 



r2X2 {(j)^{x2) + 0J(X2)) </>i(x3) " r3X3</);^(x2) (</)^(x3) - ^^(xs)) 



XE'aiQKaixU X2, X3,bub2 



(36) 



M: 



SP 
aK 



X 



327tCfM^/\/6 / dxidx2dx3 / bidbib2db2(f)B{xi,bi) 
Jo Jo 

I (X3 - 1)0;^(X2)0:^(X3) - 4r2r30j'(x2)0]^(x3) 
+r2r3X3(0j'(x2) + 0J(X2))(0^(X3) - ^^(xs)) 
+r2r3(l - X2)«(X2) - <(X2))(0^(X3) + (pKi^s)) 
XE'^{td)hna{Xi,X2, X3, 61, 62) 



+ 



X20;J^(x2)0i(x3) + X2r2r3(0j'(x2) + ^^(Xs)) (0:^(X3) - (pRi^s)) 



+r2r3(l - X3)«(X2) - <(X2))(0^(X3) + 0^(X3)) 



xKiQKaixuX2,Xs,bub2] 



(37) 
are given 



The evolution functions Ei{t) and hard functions hi appeared in Eqs. 
explicitly in Appendix B. 

If we exchange the position of K and r]^'^ in Fig. 1, we will find the corresponding decay 
amplitudes for the case of 5 — r^*^'-* transitions. Since the K and t]^'^ are all pseudoscalar 
mesons and have the similar wave functions, the decay amplitudes for new diagrams, say 

pV-A pV+A pSP fufV-A ]\/rV+A AfSP pV-A pV+A pSP ]\/rV-A ]\/rV+A J AfSP 
er) ' erj ' er} ' er) ' er) ' er) ' ar) ' ar) ' ar) ' arj ' ar) ' ciiiU. -^yJ-Q^^ 



(38) 



can be obtained from those as given in Eqs. (12811371) by the following replacements 



n(')5 



Tr) ^ Tk- 
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B 




(a) 



(b) 



(c) 



(d) 



FIG. 2. Feynman diagrams which may contribute to the B — )• Krjc decays in pQCD approach 
at leading order. 



B. B ^ Krjc at leading order 

For B — i- Krjc decays at leading order pQCD approach, the Feynman diagrams which 
may contribute are shown in Fig. O The B — )■ Kr]c decays have been studied in Ref. [23] 
in the pQCD approach, at the full leading order and with the inclusion of the partial 
NLO vertex corrections. We here recalculated theses decays and confirmed the results of 
Ref. (ill. The relevant decay amplitudes are the following 

fl /"OO 

Fri^K = SnCpM^j / dxidx3 / hidhih^dh^4)B{xi-ihi) 
Jo Jo 

:{[[(! - (1 + X3) - X3<]0^(X3) + r,il - 2x,)[^^^ixs) + ^Hx,)] 



X- 



VcK 



+r3r2j(l + 2x3)0^(x3) - (1 - 2xs)(l)lixs)]\E,{QKix,, Xs,h,h 

+2r3(l - rl)<f)^{x^)Ee{Qh',{x3, x,, h, 

32 f°° 

-^TrCpM^ / dxidx2dxs / bidbib2db2(pB{xi,bi)(j)'^^{x2) 

V6 Jo Jo 

X3(l - 2r^j0^(x3) - 2x3r3(l - O0^(x3) 
/) 



(39) 



X 



■K{t'f)h'n{^l: ^2, Xs, bi, 62) 



(40) 



where r^, = rriQ /ms, t^^ = rrir^^/mB, r^. = rric/mB- 0,^^ is the leading twist-2 part of the 
distribution amplitude for the pseudo-scalar meson r]c. The evolution function E^'\t), 
hard function hi and the scale tg, t'^ will be given in Appendix B. 

At the leading order pQCD approach, the total decay amplitude for B — )■ rjcK decay 
can then be written as: 



M{B^r],K) = F^^KU 



V*^Vcsa2 - V^^yts{a:>, -05-07 + 09 



(41) 



C. The total decays amplitudes of S — )• Krj^'^ decays 

For B^ — )• K^T] and B^ — )■ K^r] decays, by combining the contributions from all 
possible Feynman diagrams Fig. [T] and Fig. m one finds the general expressions of the 
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total decay amplitudes (here the Wilson coefficients and CKM factors are all included): 



+ {2c, + iCio) M;^-^'''Fi(0) + (^Cs - ^Cr) M^J^^'F,{4>) 

+ (^2^6 + ^Cs) M'^^'^^FM) + (a - ^^^s) Mf^'^F2(0) 

+ (^a4 - ^aio) [fBF:^^F,{cj>) + SkF^-^F.^c^) + /bFJ;-^Fi(</.)] 

+ (^Ca - ^Cg) [Mr^-^F,(0) + Ml-^F,{cl>) + M„V^Fi(</,)] 

+ [c, - \C,^ K/^F,(0) + M^+^FM) + Mr/^Fi(</,)] ] } 
+M{B^r],K)*F,{e,4>GAQ), 



(42) 



+ai/;,Fj;;-^Fi(0) + Ci [M^^^F^{<t^) + MJ;-^Fi(0) + M„V^Fi( 
-At [ ^203 - 2a5 - ^07 + ^ag^ /^^F^^ 



fs TpV—A 



/ 1 1 1 \ 

+ I as + a4 - 05 + - - 2^^io ) Jn^eK 

+ (^a, - ias) f^Ff^ + (Cs + C4 - ^Cg - ^C^o^ M^^^'^F^ 
+ (2C4 + ^Cio) M:^-^'^F,{cI>) + (c, - ^^7) M^^^'^F,{cl>) 

+ (^2^6 + ^Cs) Mf^''^Fi(0) + (^Ce - M^^'^FM 

+{a, + aio) [fsF^K^FM + f^F^'^F^icfy + /sF„V^Fi(0)] 
+(a6 + as) [/BF<fi'F2(0) + fKF^rfF,{<P) + /sFjf Fi(0)] 
+(C3 + Cg) [Mr^-^F2(0) + Mr,-^Fi(</,) + MX-^Fi(0)] 

+(a + C7) [m:/-^f2(0) + Mr/-^F,(0) + mx+-^fi(0)] ] } 

+A^(B^r/ei^)*Fe(^,0G,0Q), 



(43) 
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where A„ = V*^Vus, Ai = V^lVts, while the Wilson coefficients are defined as the combi- 
nation of Wilson coefficients Cj through the following relations: 

ai,2 = C'a,! + = + ^1^, for i = 3,5,7,9; 

a^ = Q + ^, for ^ = 4, 6, 8, 10. (44) 

The expressions of the mixing parameters (0) depend on the choice of the different 
mixing schemes: 

(i) In MS-1, i.e. the FKS r]-r]' mixing scheme, -Fi,2(0) and F[2{4>) have been defined in 
Eq. ©. 

(ii) In MS-2, i.e. the rj-rj'-G mixing scheme, Fi^2{4') and -^12(0) of ^^e form 

Fi(0) = ^p(cos + sin 6* sin 6'jAG), -^2(0) = — sin0 + sin6'cos6'jAG', 
v2 

F^(0) = ^(sin0 - cos6'sin6'jAG), ^2(0) = cos0 - cos6'cos6'iAG. (45) 

(iii) In MS-3, i.e. the third Tj-rj'-G-rjc mixing scheme, -Fi.2(0) and F[2{4') are are the 
same as given in Eq. (H5|l . For this case, the rj^ also contribute through mixing, and 
the relevant mixing parameters are 

Fc(6',0g,0q) = -sin6'sin0Gsin0Q, F^(6', 0g, 0q) = cos 6' sin0G sin0Q. (46) 

Finally, the total decay amplitudes for K^rj' and — t- K~^t]' in the pQCD 

approach at leading order can be obtained easily from Eqs. (H2l) and f H3l) by the following 
replacements 

fi -fi — s. f * 

J Tj ' J ri' 1 J q ' J Tj'i 

Fi(0) -> Fi'(0), F2(0)^F^(0), F,(0)^F,'(0). (47) 



IV. NLO CONTRIBUTIONS IN THE PQCD APPROACH 



As is well-known, the power counting rule in the pQCD factorization approach [2l[ 
is rather different from that in the QCD factorization p, 30, Hf. When compared with 
the previous LO calculations in pQCD [9], the full pQCD predictions should include the 
following NLO contributions: 



'1) We should use the Wilson coefficients Ci(mvy) at NLO level in the NDR scheme 
j28|, the NLO renormalization group (RG) evolution matrix U{t,m,a) as defined 



in Ref. [28(], and the strong coupling constant as{t) at two- loop level. 
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(2) Besides the LO hard kernel H^^\as), all the Feynman diagrams, which contribute 



to 

diagrams which contribute to H^^^ (a 
shown in Figs. 3-5, and can be classified into six types. 



in the pQCD approach, should be considered. The typical Feynman 

at the NLO level in the pQCD approach are 



(i) 



The Vertex Correction(VC): the NLO contributions from the Feynman dia 
grams as shown in Fig.3(a)-3(d), and have been evaluated ten years ago(2ll 



m 31 



m 



IV 



The Quark-Loops(QL): the NLO contributions from the quark-loops as shown 
in Fig.3(e)-3(f), the relevant analytical formulas can be found in Ref. j2l| . 

The magnetic penguins (MG): the NLO contributions from the operator Osg, 
as shown in^Fig.3(g)-3(h). These Feynman diagrams are evaluated several 



32 



years ago 

The NLO form factors(FF): i.e. the NLO contributions to the B ^ P transi- 
tion form factors F^'^^{0) with P = {K, rjq) in this paper. The typical relevant 
Feynman diagrams are shown in Fig.4, and were calculated very recently in 
Ref. fir 



The NLO contributions from the spectator diagrams as shown in Figs. 5(a)- 
5(d), which are obtained by adding a new gluon line between any two quark 
lines in Figs. 1(c)- 1(d), or by replacing the one-gluon lines with a three-gluon 
vertex in Figs. 1(c)- 1(d). Such contributions are still unknown now. 

The NLO contributions from the annihilation diagrams as shown by Fig. 5(e)- 
5(h), which are obtained by adding a new gluon line between any two quark 
lines in Figs. 1(e)- 1(h). Such contributions are also unknown now. 

The NLO contributions from the Feynman diagrams in Fig. 3, the vertex corrections, 
the quark-loops and chromo-magnetic penguins, have been evaluated several years ago 



VI 



Ref. 12 



21l . l30l |32|| . and taken into account in our previous studies for the B — )• Krj^'> decay in 



The Feynman diagrams as shown in Fig.4 can provide the NLO contributions to the 
B ^ P transition form factors and have been calculated very recently in Ref. [l3|- The 
authors of Ref.[l3| calculated the NLO corrections to the i? — t- tt transition form factors in 
the leading twist in the kx factorization theorem, and they found that the NLO part can 
provide a 20 — 30% enhancement to the LO results for the corresponding form factors. 
Since vr, K and 77*^'-* are all pseudo-scalar mesons and have the similar forms of wave 
functions, it is straight forward to extend the calculations in Ref. [13] to the cases for 
B ^ K, 1]^'^ transition form factors. In this paper, we will consider the effects of the NLO 
part of the form factors. According to general expectation, the enhanced form factors 
may lead to a larger branching ratio for B — )■ Kr]^'^ decays. 

The still missing NLO parts in the pQCD approach are the O(a^) contributions from 
non-factorizable spectator diagrams and annihilation diagrams, as illustrated by Figs. 5(a)- 
5(h). The analytical calculations for these Feynman diagrams are still absent at present. 
But it is generally believed that the NLO contributions from these Feynman diagrams 
should be small 

(i) For the non-factorizable spectator diagrams in Fig. 1(c)- 1(d), their contributions at 
leading order are strongly suppressed by the isospin symmetry and color-suppression 
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B 



K{rP) 



(a) 




(e) 



(b) 



(f) 



(c) 



(g) 



(d) 



Li/ LI LI 



(h) 



FIG. 3. The typical Feynman diagrams which provide NLO contributions to — t- Kt]^'^ decays in 
the pQCD approach: 3(a)-3(d),the vertex corrections; 3(e)-3(f), the quark-loops; and 3(g)-3(h) 
the chromo- magnetic Penguins Osg. 



with respect to the factorizable emission diagrams Fig. 1(a)- 1(b). The NLO con- 
tributions from Figs.5(a)-5(d) are higher order effects on small LO quantity, and 
therefore should be much smaller than the LO one. 

(ii) For the annihilation spectator diagrams at leading order, i.e. Figs.l(e)-l(h), they are 
power suppressed and generally much small with respect to the contribution from 
emission diagrams Fig.l(a)-l(b). The contributions from Figs.5(a)-5(d) are also the 
higher order corrections to the small quantity and therefore should be much smaller 
than its LO-part. 

In next section, we will evaluate explicitly the numerical values of the individual decay 
amplitudes corresponding to different Feynman diagrams, and will compare the size of 
every parts of the total decays amplitude for the considered decays. We try to make 
a simple and clear comparison between the contributions from different set of Feynman 
diagrams or from the different sources numerically. 



B 



\/ \/ A./ \/ 



FIG. 4. The typical Feynman diagrams which may provide NLO contributions to B ^ P form 
factors. 



A. The NLO contributions to S — )• Krj^'^ decays 



In Ref. 12[], by using the ordinary FKS r] — r]^'^ mixing scheme, we calculated the 
branching ratios and CP-violating asymmetries of the four B — )■ Kr]^'^ decays at leading 
order and partial next-to-leading order, i.e. the NLO contributions from the feynman 
diagrams in Fig. 3 were taken into account in Ref. [l2|. For the details of the calculations 
about these NLO contributions and the expressions of relevant functions, one can see 
Ref. [l3]. For the sake of reader, we here give a brief summary about those "old" NLO 
parts. 
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A. — 




FIG. 5. The typical Feynman diagrams which may provide NLO contributions to S — )• Kr/^'^ 
decays in the pQCD approach: 5(a)-5(d), the spectator diagrams; 5(e)-5(h), the annihilation 
diagrams. 



(a) Vertex corrections(VC): 

The vertex corrections to the factorizable emission diagrams, as illustrated by 
Fip. 2 (a)-2(d), have been calculated years ago in the QCD factorization appeoach[7|, 



30|,l3l|. The difference between the cases of considering or not considering the parton 



transverse momentum kr are very small and can be neglected [21|]. The NLO vertex 
corrections will be included by adding a vertex function Vi{M) to the corresponding 



Wilson coefficients aj(/x) |30l l31 



a,{fi)^a,{fi) + ^4^Cp^^^^V,{M), for z = 3,5,7,9, 
a,(/x)->a,(/x) + ^C^^^Mr,(M), for z = 4,6,8,10, (48) 

where M is the meson emitted from the weak vertex. When M is a pseudo-scalar 



meson, the vertex functions Vi{M) can be written as [2ll. |31 



V,iM) 



l2\nniL-l8 + h^J^dx(f)Ux)g{x), for ^ = 1 - 4, 9, 10, 
_l2\n^ + 6-^j;dx(j)U^)gil-x), for ^ = 5, 7, 



(49) 



6+'£lodx^Ux)h{x), 



for 



6,8, 



where /m is the decay constant of the meson M, the hard-scattering functions g{x) 
and h{x) can be found in Ref. l3 . 

(b) Quark loops(QL): 

The contribution from the so-called "quark-loops" is a kind of penguin correction 
with the four quark operators insertion, as illustrated by Fig.3(e)-3(f). We here 
include quark-loop amplitude from the operators O12 and O^^e only. The quark 
loops from Oy.io will be neglected due to their smallness. 
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For the 6 — )■ s transition, the contributions from the various quark loops are given 
by: 



H. 



q=u,c,t ql 



V2 



2n 



C^if,, I') {hp (1 - 75) T'^s) {q'YTy') , (50) 



where is the invariant mass of the gluon, which attaches the quark loops in 
Fig. 3(e)-3(f). The functions C^ifiJ^) can be found in Ref. [ij. 

The "quark-loop" contribution to the considered B — )■ Kr]^'^ decays can be written 



as 



m: 



Krj 



G p 
G 



E n^^^ Kt ^2(0) + M^^t ^1(0)1 , (51) 



q=u,c,t 



M^'? = < mnifflB >= ^ E K^^^^ ^kI m) + M^^l, F[i<t>)\ , (52) 



q=u,c,t 



where Fil{(j)) are the mixing parameters and have been defined in previous sections, 
while the decay amplitudes m]^^^ and M^^]^ are of the form 12 



M 



(9) 



Kris 



dxidx2dx3 / hidhihzdhz4>B{xi.ihi 



■ { [(1 + 0:3) 0^(X3X(X2) + (1 - 2x3) (</)^(2;3) + </)^(x3)) <P^X^2) 

+%s0x(a;3)0j(a;2) + 2ri^r^^ ((2 + X3) 0:^(a;3) - a;30x(a^3)) 0jl(a:2)] 

•^(^)(^9>^')^e(a:i,a;3,6i,63) 
+ [2rx(/.^(a;3)(/>;^Jx2) +4ri^r^.0x(a:3)0jl(a;2)] 

•^('')(t;,/")/ie(x3,xi,63,&i)}, (53) 



for B ^ K transition, and 



(9) 



76 







dxidx2dxs / bidbib^db-^ (f)B{xi,bi 



(1 + X3) 0:^(^3)0^(3:2) + (1 - 2x3) (0j;(x3) + 0J^(X3)) 0^(X2) 



+2rK0^Jx3)0^(x2) + 2r^rA' (^(2 + X3) 0^,(^3) - X30^^(X3) ) 0/r(x2) 
•E('?)(tg,/2)/,,(xi,X3,6l,63) 

+ 2r^0j,(a:3)0K(a:2) +4r^ri^0j'^(x3)0^(x2) 

•£;('')(t;,r2)/i,(x3,xi, 63,^1)}, 



(54) 



for B ^ 7] transition. Here = ml/mB and r^^ = rriQ/mB- The expressions 
of evolution functions E^'^\tg,P) and E^'^^t'^J''^), as well as the hard functions 



/ie(a;i, X3, 61, 63) and /;,e(a;3, Xi, 63, 61) can be found in Ref. [12 
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(c) Chromo-magnetic pengmns(MP): 

This is another kind of penguin correction but with the Osg insertion. The corre- 



sponding effective weak Hamiltonian for b — > sg transition is of the form [21 

Gf 9. 



TTcmp 



V2 8tt' 



(55) 



where the effective Wilson coefficient Cl{^ = Csn + 21 



The total chromo-magnetic penguin contribution to the considered B — )■ Ki]^'^ de- 
cays can be written as 



M 



(cmp) 
Kri 



j^Acmp) 



< K7]'\nijf\B > 



V2 
Gp 

71 



A, 



(9) 



(56) 
(57) 



where the mixing parameters F^\{(1)) have been defined in Sec. II. The decay ampli- 
tudes and M^^]^ are obtained by evaluating the Feynman diagrams Fig. 3(g)- 
3(h)[l2|: 



(9) 



1 noo 

dxidx2 dx'i 
Jo 



bidbib2db2b3db3 0b(xi, bi) 

Jo Jo 

•||2(-l+X3)0;^^(x2)0/,(x3)+r,^X2(l + X3)[-30j;(x2)+0j,(x2)]0i(x3) 

+rK[{-S + 2x3 + xDcP'^ixs) + (-1 + 2a;3 - xI)(PI{x3)](PI{x2) 

+3r^,ri^[(-l - a;2 + 2:3 + 2x2X3)0^(^3) + (1 - X2 - X3 + 2x2X3)0|^(x3)]0j'^(x2) 



+Vk[{-1 + X2 + X3 - 2X2X3)0^(X3) + (1 + X2 - X3 - 2x2X3)0^(x3)]0j'^ (X2) 

■Eg{tg)hg{A, B, G, ^l, 63, 62, ^ 3) 

-[4^rK(l)^^{x2)(j)K{x3) + 2r^^ri^X20^(x3)(3(/)^J(x2) - 0j(x2))] 

■Eg{t'g)hg{A', B', G', 63, bi, 62, xi)}, (58) 



M 



(g) 



1 noo 

dxidx2 dxs 
Jo 



bidbib2db2b3db3 0b(xi, bi 



■{{ 



2(-l + X3)0^(x2)04(a^3) + rxX2(l + X3)[-30^(X2) + 0^(X2)]04 (^s) 



+r^[(-3 + 2x3 + a;2)0^^(x3) + (-1 + 2x3 - a;2)0j;(x3)]0i(x2) 



+3r^ri^[(-l - X2 + X3 + 2x2X3)0^^ (X3) + (1 - X2 - X3 + 2x2X3)0^^ (x3)]0^(x2) 
+rnrK[{-l + X2 + X3 - 2x2X3)0^^ (X3) + (1 + X2 - X3 - 2x2X3)0^^ (x3)]0^(x2)| 

■Eg{tg)hg{A,B,G, b^MM^X^) 

-[4r^0;^(x2)0j' (X3) + 2r^ri^X20j' (X3)(30^(x2) - 0^(a;2))] 



■Eg{t'g)hg{A' ,B' ,G' MMM.Xi)].. 



(59) 
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where = ml/niB, 



evolution functions Eg and hard functions hg in Eqs. 



Ref. [ii 



mf^/mB, tk = ^0 /'nT'B- The exphcit expressions of the 

can be found easily in 



B. The form factors at NLO level 



As mentioned in introduction, Li et al. derived the /cT-dependent NLO hard kernel 
H^^^ for the 5 — t- tt transition form factor IJ]. We here extend their results for B ir 
form factors to the ones for B ^ K and B — )• {Tjq^rjs) transitions, under the assumption 
of S'f/(3) flavor symmetry. As given in Eq.(56) of Ref. [13], the NLO hard kernel H^^^ can 
be written as 



47r 



21 



m 



B 



m 



— In — ( In — ^ + — ^ 1 In —5- + — ln^(xia;3) + - In^ x 



13 



1 

'4 



m 



+- In Xi In X3 + ( 2 In + —Inr] — - ) In Xi + ( - In 77 — ) In X3 



7 



7 



m 



B 



16 
7 



+ 



15 7 



In r? In r? In 

16 " ' " 



m 



B 



3 In 



m 



B 



+ 2 + 



101 



219 

vr" H 

48 16 



(60) 



where the scale Ci = 25777,^, rj = 1 — {pi — PzY /rn^ is the energy fraction carried by the 
meson which picks up the spectator quark. For B — )■ Krjq, Krjs decays, the large recoil 
region corresponds to the energy fraction r] ~ 0(1). For B — )■ Krjc-, ti{K) ~ (1 — r^^). iXf 
is the factorization scale, which is set to be the hard scales 



max 



(ya^mB, 1/61, 1/63), or = max(yxlr/ m^, I/61, 1/&3), 



(61) 



corresponding to the largest ene rgy scales in Fig. [T](a) and[T](b), respectively. The renor- 
malization scale /i is defined as 14 



= ts(yUf) = <^ Exp 



ci + (ln:^ + ^ lln-^ 



m 



B J 



with the coefficients 



1 -^2 



2/21 



(62) 



Cl 



C2 



^--^lnr/)lnr/ + ^ln^(31n^ + 2) - ^tt" - 



16 



m 



m 



B 



101 , 219 



48 



16 



2 In 



7 



mi, . , 

3 

+ 5- 



C3 = --ln?7 

o 

At the NLO level, the hard kernel function H can then be written as 
H = H^'\a,) + H^^\ai) = [I + F{x,, X3, /x/, v, Ci)] H^^^ 



(63) 
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C. The NLO contributions to S — )• Krjc decays 



For B — 7- rjcK decays, the NLO contributions include two parts: (a) the NLO vertex 
corrections to these decays, which has been taken into account in Ref. [2^; and (b) the 
NLO contributions to the B K transition form factors, which is newly known NLO 
part. 

Since the emitted meson is rjc = cc, the soft and collinear infrared divergences of the 
four vertex correction diagrams will be canceled each other. So these vertex corrections 
can be calculated without considering the transverse momentum effects of the quark at 
the end-point region, the same way as in the collinear factorrization theorem. 

The NLO vertex corrections can be included through the re-definition of the Wilson 
coefficients: 

a2^a2 + ^Cp^{-18 + 12 In ^ + /,), 

a,^a, + ^C^%i(-18 + 121n^ + /,), for z = 3,9, 
Air 6 jj, 



%^a,-^C^^^^iF^(-6 + 121n^ + /,), /or j = 5, 7, (64) 



^C^^kiM(_6 + 121n!^^ 
47r 3 n 

with the function // is of the form 

^ 2^2Nl f\ . ,r3(l-2x), ,^ , . , 2^(1 -x)] 

// = / dx(f)l (x) In X + 3 In 1 - z) - itt) + — ^ , 65 

U Jo I l-x l-zx \ 

where z = m^^/m\. 

The NLO contributions to the B ^ K transition form factors can be included for the 
B — 7- Krjc decay in the same way for B — )■ Kr]^'^ decays. 



V. NUMERICAL RESULTS AND DISCUSSIONS 
A. Input parameters 



We use the following input parameters (lOl . in the numerical calculations(all masses 
and decay constants in units of GeV) 

/b = 0.21, fK = 0.16, fr,^ = 0.4874, = 0.5475, m^/ = 0.9578, 

TTiKO = 0.498, mK+ = 0.494GeV, moK = 1-7, Mb = 5.28, 

rrih = 4.8, = 1.5, m^^ = 2.98, 

Mw = 80.41, me = 1.376, r^o = 1.53ps, tb+ = 1.638ps. (66) 

For the CKM quark-mixing matrix, we adopt the Wolfenstein parametrization as given 
in Ref. M,^- 



= X = 0.2246, \Vub\ = 3.61 x 10~^ 
= 0.9748, Vcb = 0.04197, 

Vts = -0.042, Vtb ^ 1.0, (67) 

with the CKM parameters: A = 0.2246±0.0011, A = 0.832±0.017, p = 0.130±0.018, f] = 
0.350 ±0.013. 





0.9748, 


Vcd = 


-0.225, 


Vtd\ = 


8.8 X 10 
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B. Form factors at LO and NLO level 



We firstly calculate and present the pQCD predictions for the form factors at zero 
momentum transfer for B — t- Ki]^'^ decays at LO and NLO level respectively. In the 
calculation, we consider three different mixing schemes respectively. 

In this paper the form factors F^^^{0) and Fq^^ (0) are defined as 

Fo^^''(O) = cos0Fo^^^'(O),, 

F,^^'\0)=sm{cf>)F,^^"^{0)j, (68) 

in the ordinary FKS t] — rj' mixing scheme, and 

F(f ^"(0) = [cos (j) + sme sin ^^(1 - cos (^g)] F(f ^"'(0)//, 
Fq^"'{0) = [sin (f)- cose sme,{l- cos (f)G)]F^^'^'{0)jj, (69) 

in the MS-2 and MS-3 mixing schemes. One should note that the form factor Fq'^^''(0)j is 
different from Fq^^''{0)jj since some relevant parameters in the distribution amplitudes, 
such as Prf^ = 2mq/mqq, are different in different mixing schemes. The pQCD predictions 
for the numerical values of the form factors for three different mixing schemes are all listed 
in Table I, and they are obtained by using the central values of all input parameters. For 
the relevant mixing angles, we take = 39.3° in MS-1 mixing scheme; while we take 
Oi = 54.7°, 9 = -11°, (f) = e + e^ = 43.7° and 0g = 12° in both MS-2 and MS-3 mixing 
scheme. The error comes from the uncertainty of Ub = 0.40 ± 0.04 GeV. 



TABLE L The LO and NLO pQCD predictions for the form factors of the B — t- Kt]^'^ decays 
for three different mixing schemes. 



Form factors 


MS 


LO NLO 




1 


0.20 0.26 ±0.04 




2,3 


0.28 0.33 ±0.06 


Fo^-^(O) 


all 


0.37 OAStlf^ 



From the numerical values of the form factors in Table. HI we can see that (a) the form 
factors are the same for the rj-rj'-G mixing scheme and the rj-ri'-G-ric mixing scheme, since 
the r]c component in the rj-rj'-G-rjc mixing scheme does not affect the evaluation of the 
form factors Fq^^ (0) and F^^^(O); and (b) the NLO contributions also provide ~ 20% 
enhancements to the corresponding form factors. 

C. Br{B — )■ Kf]^'^) in r]-r]' mixing scheme 

In the B-rest frame, the branching ratio of a general B — )• M2M3 decay can be written 

as 

Br{B ^ M2M3) = TB — ^ X \MiB ^ MaMs)!' , (70) 
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where tb is the hfetime of the B meson, x f« 1 is the phase space factor and equals to 
unit when the masses of final state light mesons are neglected. 

Using the input parameters and the wave functions as given in previous sections, it 
is straightforward to calculate the CP-averaged branching ratios for the considered four 
B — )■ Kt]^'^ decays in different mixing schemes. For the case of the ordinary r] — r]' mixing 
scheme, the pQCD predictions are listed in Table. Ull where the label "NLO-l" refers 
to the pQCD predictions with the inclusion of the same set of NLO contributions as in 



Ref. [12|. The label NLO in Table III means all currently known NLO contributions 
are included, especially the NLO part of the form factor obtained by evaluating the 
Feynman diagrams as shown in Fig. 4 For comparison, we also list the corresponding 
experimental measurements fl3] and the theoretical predictions in the pQCD approach 
[3 and in the QCDF approach (3l|. 



TABLE IL The pQCD predictions for the branching ratios (in unit of 10~^) in the ordinary ry-?]' 
mixing scheme with (j) = 39.3°. The label NLO means all currently known NLO contributions 
are included. 



Channel 


LO NLO-l NLO 


DataflO] 


pQCD[12] QCDF[31] 


B+ K+7] 
B+ K+7]' 


2.12 2.76 2.62if:^ 
27.9 48.3 57.2+23-7 
3.83 3.78 4.0l|^ 
30.3 49.8 58.7tlti 


1 r,q+0.27 
i.ZO_g 24 

66.1 ±3.1 

9 /I +0.22 
^•^-0.21 

71.1 ±2.6 


9 1+2.6 1 1+2.4 
^■-■--l.S ^-^-1.5 

^'^•'^-10.6 ^0-^-22.0 
q 9+3.2 1 Q+3.0 
^-^-1.9 
c;i n+18 40 1+45.2 
^•^•0-10.9 49.i_23.6 



Of course, the NLO pQCD predictions for the CP-averaged branching ratios still have 
large theoretical uncertainties. If we take into account the effects of the uncertainties of 
the main input parameters, we find that 

Br{B° ^ K^T]) 
Br^B"" ^ K%') 
Br{ B+ K+T]) 
Br{ fi+ ^ K+T]') 

where the major errors are induced by the uncertainties of uji, = 0.4 ± 0.04 GeV, = 
0.13 ± 0.03 GeV, fs = 0.21 ± 0.02GeV and Gegenbauer moment = 0.44 ± 0.22 ( 
here denotes a^'' or al"), respectively. The total theoretical errors in the NLO pQCD 
predictions as shown in the fourth column of Table II are obtained by adding the four 
individual theoretical errors in quadrature, from the numerical results as given in Eq. (!7Ti) 
and Table II we fnd the following points 

(i) By comparing the predictions as listed in "NLO-l" column and "NLO" column, one 
can see that the inclusion of the NLO contributions to the form factors can provide 
about 18% enhancement to Br{B — )■ Kt]'), the gap between the pQCD predictions 
and the measured values therefore become narrow effectively, but there is still a 
small difference between the central values of the pQCD predictions and the data. 



= [2.62^^:?^(u;,)l?:^^(^^)+o.52(^^)+i.37(^.)] ^ 10-6^ 

= [57.2ll?iK)i^2g^K)llJ:l(/^)l^:l^(a^)] x lo-^ 

= [3.97^i:?3^K)^?|^K)^°:^^(/^)l}i^(al)] X lo-^ 

= [58.7^i?:^(u;,)1^^2SK)liJ:e^(/^)l|}^(a^)] x 10-^ (71) 
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(ii) For B — )■ Kr] decays, however, the pQCD predictions for their branching ratios 
remain basically unchanged after the inclusion of NLO part of the form factors. 
Although the NLO pQCD predictions for Br{B — )■ Kr]) are consistent with the data 
within one standard deviation, but the central values of the NLO pQCD predictions 
are still larger than the measured values by almost a factor of 2. 

(iii) The pQCD predictions as given in the NLO-1 column agree well with those presented 
in Ref. [l2[(i.e. the results as listed in the "pQCD" column of Table II), the small 
differences are induced by the variations of some input parameters, such as the 
CKM matrix elements. 

(iv) Although the NLO pQCD predictions for Br{B Kif) and Br{B — ?■ Krj') are 
consistent with the data within one standard deviation, but we here can not provide 
a good interpretation for the observed pattern of Br{B — )■ kr]^'^) in the FKS r]-r]' 
mixing scheme. 



D. Br{B Kr]'^'^) in r]-r]'-G mixing scheme 

In the rj-ri'-G mixing scheme, by using the input parameters and the wave functions 
as given in previous sections and fixing the mixing parameters 6 = — 11°,0 = 43.7° and 
(pG = 12°, we find the LO and NLO pQCD predictions for the CP-averaged branching 
ratios as listed in Table III. As a comparison, we also show the measured values and the 
QCDF predictions as given in Ref. |3l| in last two columns of the Table III. 

TABLE III. The same as in Table I but for the case of the r]-r]'-G mixing scheme with the 
mixing parameters 4> = 43.7°, 6 = -11°, 61^ = 54.7° and (pc = 12°. 



Channel 


LO 


NLO-1 NLO 


Data[10l 


QCDF[31] 


B^- 




0.90 


1.15 


i.i3t?:° 


-I 00+0.27 
l.ZO_Q 24 


-L-J^-l.S 


5°- 


-> K^T]' 


35.2 


57.4 


66.51?^:^ 


66.1 ±3.1 


46 5+41-9 


B+ - 


K+rj 


1.98 


2.10 


2.36+f^ 


Q Qc+0.22 
Z.OD_o.21 


1 q+3.0 


B+- 


K+r]' 


38.9 


58.3 


67.3t?^.o 


71.1 ±2.6 


AO 1+45.2 
^^■-'--23.6 



The NLO pQCD predictions with the inclusion of the major theoretical errors are the 
following 

Bri B' ^ K%) = [ll3l°j°(a;,)lJ|?K)+°:i(/^)l°|^(a^)] x 10-^ 
Br{ ^ ir%') = [66.5lJ|^(a;,)l|^K)l}|?(/^)l|?(a^)] x 10-^ 
Bri ^ K^r^) = [2MtrJiiu^,)tlf,im,n^^^ x 10-^ 

Bri ^ K^v') = [67.3li|^(a;,)l^°4'K)l}|^(/^)l?i(a^)] x 10"^ (72) 

Analogous to the case of MS-1, the major theoretical errors in MS-2 mixing scheme 
are still induced by the uncertainties of the input parameters: Ub = 0.4 ± 0.04 GeV, 
= 0.11 ± 0.02 GeV, fe = 0.21 ± 0.02GeV and Gegenbauer moment = 0.44 ± 0.22. 
The total theoretical errors of the NLO pQCD predictions in the fourth column of Table 
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Ill are obtained by adding the four individual theoretical errors in quadrature. One can 
see from the numerical results in Eq. f l72]) and Table III that 

(i) In the rj-ri'-G mixing scheme, the glueball part play an important rule in improving 
the agreement between the pQCD predictions with the data. Even at the leading 
order, the pQCD predictions for Br{B Kt]') {Br{B — )■ Kr])) become larger 
(smaller) than those in the case of MS-1, the corresponding changes are what we 
want to interpret the data. 

(ii) For 5° K^r]' {B+ K+t]') decay, the NLO contribution provide a 89% (73%) 
enhancement to its branching ratio with respect to the LO prediction. The NLO part 
of the form factors along provide a 23% enhancement to the corresponding branching 
ratios. The NLO pQCD predictions for both Br{B^ K^r]') and Br{B+ K+r]') 
are now in full agreement with the data. 

(iii) For both B^ — t- K^t] and B^ — t- K^t] decays, the NLO enhancements are small 
in size, the agreement between the pQCD predictions and the data also improved 
effectively due to the inclusion of the all known NLO contributions. 

(iv) For all four B — )■ Kri'^''> decays in consideration, the differences between the numer- 
ical values as listed in the NLO-1 column and NLO column in Table III show the 
effects of the inclusion of the NLO part of the form factors. It is easy to see that 
the NLO pQCD predictions for the branching ratios become in perfect agreement 
with the measured values due to the contribution from the guleball component in 
the rj-rj'-G mixing scheme and the inclusion of the NLO contributions. 

E. Br{B Kr]^'^) in ^^r]-r]'-G-r]c" mixing scheme 

In the ^^rj-ri'-G-rjc" mixing scheme, by using the input parameters and the wave func- 
tions as given in previous sections and fixing the mixing parameters 6 = —11°, (p = 43.7°, 
(pG = 12° and (pq = 11°, we find the LO and NLO pQCD predictions for the CP-averaged 
branching ratios as listed in Table IV. 

TABLE IV. The same as in Table I but for the case of the "t/-7/'-G-?7c" mixing scheme. The 
label NLO means all known NLO contributions are included. 



Channel 


LO NLO-1 NLO 


Data[10] QCDF[31j 


B° K^T]' 
B+ K+r] 
B+ -+ K+rj' 


0.67 0.87 0.82l:J| 
43.5 55.6 64.8^^^;^ 
1.50 2.00 2.19tli 
51.7 56.2 65.21^7-0 


1 r,o+0.27 1 -1+2.4 
i.ZO_Q24 -L.-L-LS 

66.1 ±3.1 46.5^^^;^ 

r, qC + 0.22 -1 Q + 3.0 

z.oo_Q2i -L-y-i.g 

71.1 ±2.6 



In the ^^rj-rj'-G-ric" mixing scheme, the contributions from the decay chain B — )■ Kr/c — )■ 
Kt]^'") will be included. The NLO pQCD predictions for the CP-averaged branching ratios 
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with the major theoretical errors are of the form 

Bri ^ = [0.82l°:?^(c.,)lJ:^^K)l°:JI(/«)l°:i(a^)] X lo-^ 
Br{ 5° ^ KY) = [64.8l?l:?(u;,)l^°53K)li?:^(/^)] x 10~', 
Bri ^ K+v) = [2.19l°:^^(a;,)lg:°^K)l°;l^(/^)lJi^(a^)] x 10-^ 
i?r( 5+ ^ ir+V) = [65.2^?^:?K)l^°/K)+i?:°(/^)] X 10-^ (73) 

Analogous to the cases of r]-T]' and rj-rj'-G mixing schemes, the major errors here are also 
induced by the uncertainties of cUb, m^, fs and Gegenbauer coefficient al, respectively. For 
B — 7- Kt]' decays, however, the error induced by the uncertainty of al = 0.44±0.22 are very 
small and have been neglected. The total theoretical errors of the NLO pQCD predictions 
in the fourth column of Table IV are obtained by adding the individual theoretical errors 
in quadrature. 

From Table IV, one can see that the NLO pQCD predictions for Br{B Ki]') also 
agree well with the data. For B — i- Krj decays, although the central values of the NLO 
pQCD predictions for Br{B — >■ Kr]) are a little smaller than the measured values, but 
they are still consistent with the data within one standard deviation. Since the values of 
the relevant mixing parameters Fc{9, (pa, 4>q) and Flj{9, (pc, <t>Q) as defined in Eq. (1471) are 
all very small, 

Fc{0, <Pg, <Pq) = 0.0076, F'c{e, 0g, <Pq) = 0-039, (74) 

for (6*, 0G, 0q) = (— 11°7 12°, 11°), the rjc contributions to the B — )■ Krf^'^ decays are indeed 
very small. 



F. CP-violating asymmetries 

Now we turn to the evaluations of the CP-violating asymmetries of i? — Krj^'^ decays 
in pQCD approach. For B^ — )■ K^i]^'^ decays, the direct CP- violating asymmetries Acp 
can be defined as: 

... ^ ns' ^ /) - nB^ ^ /) ^ \Mj?-\Mf\- 



Using the input parameters and the wave functions as given in previous sections, it is 
easy to calculate the direct CP-violating asjTiimetries for the considered decays, which 
are listed in Table El The major theoretical errors as given in Table |V] are induced by 
the uncertainties of input parameters of Ub, irts and a^. As a comparis on, we also list 



currently available data [lO| and the corresponding QCDF predictions |3l|. The label 
NLO means that all known NLO contributions are taken into account. For B^ — )■ K^rj 
decays, there is a large direct CP asymmetry {A^'^p ), due to the destructive interference 
between the penguin amplitude and the tree amplitude. 

From the pQCD predictions and the relevant data as listed in Table V, one can see 
that 

(i) For B^ — )■ K^T] decays, the LO pQCD predictions for A'^}, in all three mixing 
schemes have an opposite sign with the measured value. The inclusion of the NLO 
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TABLE V. The LO and NLO pQCD predictions for the direct CP asymmetries A^piB^ 
K^rj) and A'qp{B^ — t- K^rj') in the three different mixing schemes (in units of 10~^). 



Mode 



MS 



LO NLO 



Data[10] QCDFfSl] 



10.0 
31.1 
42.4 



-25.2l?-^(c^,)l}0:|(m.)+8-0 



-22.9l^i(^^)+ii.6(^^) 
-2.8"'"" 



-5.2 
-8.5 
-8.5 



16 

;-H.9 

'-3.3 



12.2 
6.3 
7.6 
8.4/ »? 
3.3V"'2 



{4) 



-37 ±8 -18.9 



-29.0 
-30.0 



-10.4 
-12.2 
-9.0 



-4.4^ 



0.6 
-0.8 



0.8 
-1-0.9/ 



'-1.3 
1+1.5 
'-1.5 



(4) 



1 3+^-6 



_q n+10.6 
^•U-16.2 



-2.3 



+1.1 
1.1 



+0.4, 
0.4' 



contributions changed the sign of the pQCD prediction for A^p, the NLO pQCD 
predictions for Af^plB^ — )■ K^rj) in the cases of MS-1 and MS-2 are now becoming 
consistent with the data within one standard deviation. In the ^^rj-'q'-G" mixing 
scheme, for example, the NLO pQCD predictions are of the form 

At^UB^ ^ K^V) = -22.9^}^:? X 10-^ 

^ K^V') = -5.5^};^ X 10-^ (76) 

where the individual errors as shown in Table V are added in quadrature. 

In the case of MS-3, however, the pQCD prediction for Aqp{B^ — )■ K^rj) changed 
its sign after the inclusion of the NLO contributions. The NLO pQCD prediction 
is of the form 

A^SpiB^ ^ K^V') = -2.81^°/ X 10-2, (77) 

which is still much smaller than the measured value in magnitude. There is a clear 
difference between the pQCD prediction and the data for Aqp{B^ — > K^rj') in the 
^^rj-'q'-G-'qc" mixing scheme. 



IS 



(iii) For K^i]' decays, the measured value of Acp{K^r]') = l.Stl'^ ^ 10"^ 

consistent with zero. The NLO pQCD predictions in three different mixing schemes 
agree well with the data within one standard deviation, while the consistency be- 
tween the pQCD predictions and the data are improved effectively by the inclusion 
of the NLO contributions. 

As to the CP- violating asymmetries for the neutral decays B^ K'^r]^'^ , the effects of 
B^ — 5° mixing should be considered. The CP- violating asymmetry of B^{B^) K'^r]^'^ 
decays are time dependent and can be defined as 

r^o iAt) -Foo ..(At) 
AcP ^ ' ' = Cf cos(AmAt) + Sj sin(AmAt), (78) 
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where Am is the mass difference between the two mass eigenstates, At = tcp — Uag 

is the time difference between the tagged 5° {B ) and the accompanying B (5°) with 
opposite b flavor decaying to the final CP-eigenstate fcp at the time tcp- The direct and 
mixing induced CP- violating asymmetries Cf ( or in term of Belle Collaboration) and 
Sf can be written as 

Adir n _ l-^l ~ 1 litnix q _ 2/m(A) , , 

■^CP - W = ^ ^ U|2' ^CP - = ^ ^ U|2' 



with the CP- violating parameter A 



By integrating the time variable t, one finds the total CP asymmetries for B^ — t- K'^t]^'^ 
decays, 

/[tot Adir 1 A mix 



where x = Am/T = 0.774 H 



In Table IVIt we show the LO and NLO pQCD predictions for the direct, the mixing- 
induced, and the total CP asymmetries for B^ — t- K'^t]^'^ decays in three different mixing 
schemes. Analogous to the NLO pQCD predictions for the branching ratios, the label 
"NLO" here means the inclusion of all currently known NLO contributions are taken into 
account. 

From the pQCD predictions and currently available experimental measurements for 
the CP-violating asymmetries of B^ — )■ K'^i]^'^ decays, one can see that: 

(i) Unlike the cases for the branching ratios, the pQCD predictions for the CP-violating 
asymmetries of the neutral B'^ — )■ K^t]^'^ decays are not sensitive to both the NLO 
contributions and the choice of the mixing schemes. 

(ii) The NLO pQCD predictions for ^ iTgr/') and ^ iT^r/') have 
small theoretical error and agree very well with the measured values: 

A'^SUB'' ^ K^sV') = (3.3 ± 0.3(theory)) x 10-^ 

A^'^{B° K^r]') = (70.3 ± 0.5(theory)) x lO'^, (82) 

while the measured values are (1 — 9)% and (64 ± 11)%, respectively. 

(iii) The pQCD predictions of A'^UB^ ^ ^s^) ~ -16% and A'S'piB^ ^ K'^sV) ~ 67% 
will be tested by the LHCb and the forthcoming Super-B experiments. 
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TABLE VI. The pQCD predictions for the direct-, mixing- induced and total CP asymmetries 
(in units of 10"'^) for — )• K^rj^'^ decays in three different mixing schemes, and the world 
averages as given by HFAG [lo| . 



Mode 


MS 


LO 


NLO 


Data 






1 


-4.6 


-ll.l+°-?(wfe)+2-9(m,)+2-6(a^) 

— U.7 V ^ / — z.U \ ^ y —o.o\ 2 / 






^ ^S^l 


2 


-6.6 


-16.0+°:^(a;,)+^|(m,)l^i3^(a^) 


- 






3 


-7.8 


-19.4+°-^(wfe)+?fi^(m,)+Ifi^(a^) 








1 


69.3 


66.3+^-^(a;fe)+^1(m,)+^-2(a^) 




/[mix ( rjO 




2 


69.9 


66.7+j-^(a;b)+2-6(m,)+2-2(a^) 

— 1.5 V '''' — D.5 V / — 5.2 V 2 / 








3 


70.3 


69.5+2-i(a;fe)+°i(m,)+2-6(a^) 








1 


30.6 


25.2+0-3(a;b)+3-0(„iJ+2.7(o^) 

— U.2V "^z — 2.yv / — 4.UV 2/ 




A^^%{B^ - 

•^CP^ 




2 


29.7 


22.3+°-2(tJ;,)+4-8( )+4.4(^) 








3 


29.2 


21.6il-^K)t5f^(^^)+4.9(„^) 









1 


1.1 


3.4l0:2(u;,)10i(m,)l°:l(a^) 




A'^'piB'^ - 




2 
3 


1.0 
0.9 


3.3l0i(u;,)10:J(m,)10:J(a^) 


1±9 






1 


70.7 


69.8^°:l(..,)^0:?(m,)+°:2(a^) 






-> Kli) 


2 


70.8 


70.0t°:i(c.,)tO:}(m,)l°:J(a^) 


64± 11 






3 


70.8 


70.5t°:}K)tO:J(m,)t°:J(al) 








1 


34.9 


35.9t0:2(..,)t0:0(m,)t0:0(a^) 




A'gUB'^ - 




2 


34.9 


36.0t0i(a;,)tH("^.)lS:S(al) 








3 


34.8 


36.3i°:l(a;,)i0:l(m,)l°:l(al) 





G. Relative strength of the contributions from different sources 

In the pQCD approach at leading order, we generally expect that (a) the factorizable 
emission diagrams Figs. 1(a) and 1(b) provide the dominant contribution to the considered 
B Kr]^'^; (b) the non-factorizable spectator diagrams Figs. 1(c) and 1(d) are strongly 
suppressed by both the isospin cancelation and the color-suppression and therefore play a 
minor role; and (c) the annihilation diagrams Fig.l(e)-l(h) are generally power suppressed 
in magnitude, but may provide a large strong phase to produce large CP-violating asym- 
metries for some decay modes. 

In the pQCD approach at next-to-leading order, as presented in previous sections, we 
assume that those currently known NLO contributions are the dominant part of the NLO 
contribution, while the still missing part from the spectator and annihilation diagrams 
Figs.5(a)-5(h) is small in size and can be neglected safely. Of course, this assumption 
should be examined properly before the calculations for those missed part become avail- 
able. For this purpose, we try to check about the relative strength of the LO or NLO 
contributions coming from different sources. 

If the LO contributions ( oc 0(a<j)) from Figs.l(c)-l(d) and Figs.l(e)-l(h) are already 
much smaller in size when compared with those from the emission diagrams Figs. 1(a)- 
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1(b), it is reasonable for us to assume that the still missing next-to- leading order 0{a1) 
corrections coming from Figs.5(a)-5(h) should be smaller than their counterparts at lead- 
ing order, and therefore much smaller than those dominant LO contribution and therefore 
can be neglected safely. 

In order to check these general expectations or assumptions correct or not, we here 
will firstly decompose the LO decay amplitude Ai lo into different parts according to the 
corresponding Feynman diagrams, and then make numerical evaluations for each part 
and compare their magnitude directly. We try to make a simple and clear numerical 
comparison for the contributions from different sources 

For — > K~^ri decay in the rj-rj' mixing scheme, for example, its decay amplitude 
A4.{B~^ — )■ K^rj) at leading order as given in Eq. gSDEl can be written as a sum of three 
parts 

Mlo{B+ K^V) = M^^^K+ri) + M'+\K+7]) + TW"""^^'^^), (83) 

where the decay amplitude A^"+^ are obtained by evaluating the dominant emission di- 
agrams Figs. 1(a)- 1(b), }A^^'^ refers to the LO contribution from the spectator diagrams 
Figs.l(c)-l(d), while denotes the LO contribution from the annihilation diagrams 

Figs.l(e)-l(h). 

By using the central values of the input parameters and the relevant wave functions, 
we make the numerical calculations and find the results (in unit of 10~^) 

Mlo{B+ K^V) = ,-1-76 - i0.37 + 0.065 - i0.14 + 0.01 + iO.57 

= -1.67-i0.062, (84) 

it is easy to see that (a) A^'*"'"* is indeed large and dominant; (b) Ai'^'^'^ is negligibly 
smaller than both A^"+'' and and (c) the real part of ^Vf^"™ is close to zero, but 

its imaginary is large and interferes destructively with A^""'"''. 

Since the branching ratio of the considered decays are proportional to the square of 
the decays amplitude |A^P, as shown by Eq. (ITOl) . we can define the relative strength of 
the individual contribution from different sources as the ratio Rlo'- 

RLo{K+r]) = iTW^+^p : \M^+d.\2 . ^^anm^2 . \Mlo\^ = 3.23 : 0.02 : 0.33 : 2.79. (85) 

One can see directly from above numbers that the contributions from Ai'^'^'^ is less than 1% 
and can be neglected safely, while the contributions from _A/1'*""* is also small in magnitude, 
around 10% of the dominant contribution from emission diagram Figs.l(a)-l(b). 

Using the same methods, we make the similar decompositions and numerical calcu- 
lations for the remaining three decay modes 5° — )■ K'^t]^'^ and — )■ K~^r]', and find 
the numerical values of the decay amplitudes and the relative strength. We make the 
calculations in both the rj-i]' and rj-rj' -G mixing scheme and show all the numerical results 
in Table IVIIi For the case of each mixing scheme we use the same input parameters as 
those used in the calculation for branching ratios in Sec.V-C and V-D , respectively. 

From the numerical results as shown in Table IVIII we find the following points: 



^ In the 77-77' and rj-rj' -G mixing scheme, the last term M.{B VcK) ■ Fc{9, (pC 0q) in Eq. (|43|) is absent. 
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TABLE VII. The LO pQCD predictions for the numerical values (in unit of 10~^) of the indi- 
vidual and total decay amplitudes of — t- + r/^'^ and — t- K^rj^'^ decays, and in the 77-77' 
and r]-r]'-G mixing scheme. 



Decay 


MS 








Mlo 


Rlo 




1 

2 


-1.30 - iOM 
-0.80 + iom 


0.06-i0.11 
0.03-i0.07 


-0.06-FiO.53 
0.01 -FiO.54 


-1.30 + i0.47 
-0.76-hi0.51 


1.69 : 0.01 : 0.29 : 1.91 
0.63 : 0.01 : 0.30 : 0.83 


K^i 


1 

2 


3.42 - ii0.03 
4.12 + i0.03 


0.25 - iO.47 
0.24 - iO.45 


-0.07-i2.85 
-0.03 -i2. 94 


3.40 - i3.29 
4.32 - i3.37 


11.7 : 0.29 : 8.1 : 22.4 
17.0 : 0.27 : 8.6 : 30.0 


K+T] 


1 

2 


-l.76-i0.37 
-1.18 - fO.52 


0.07-i0.14 
0.03 - iom 


0.03-Fi0.57 
O.lO-FzO.58 


-1.67-hi0.06 
-1.05 - i0.03 


3.23 : 0.02 : 0.33 : 2.79 
1.66 : 0.01 : 0.35 : 1.10 


K+T]' 


1 

2 


3.65 - i0.30 
4.44 - iO.49 


0.26 - iOM 
0.25 -ii0.05 


-0.41 -i2. 99 
-0.38-i3.07 


3.50-i3.83 
4.31 -i4.08 


13.4 : 0.36 : 9.1 : 26.9 
20.0 : 0.33 : 9.6 : 32.2 



(i) For all the four considered decays, the factorizable emission diagrams Figs.l(a)-l(b) 
provide the dominant contribution to the branching ratios. In both MS-1 and MS-2 
mixing schemes, we have 



\M 



c+d\2 



\M 



a+b\2 



< 0.01, 



I anni 1 2 



(86) 



for the two B Kr] decays, and 
for the two B — )■ Krj' decays. 



< 0.03, 



\M 



anni 1 2 



\M 



a+b\2 



< 0.7, 



^7) 



(ii) For all the four considered decays, one can see from Eqs. (186|87p that the contribution 
from the spectator diagrams Figs.l(c)-l(d) is very small in size. 



\M 



c+d\2 



\M 



LO\2 



< 0.03, 



and therefore can be neglected safely, which is consistent with the general expec- 
tation. Since the LO part Ai^^"^ is already negligibly small, it is reasonable for 
us to neglect the corresponding higher order NLO contribution Ai 
corresponding spectator diagrams Figs.5(a)-5(d). 



from the 



(iii) For all the four considered decays, the real parts of ^A/f"""* are very small, but their 
imaginary parts are relatively large, this leads to a large strong phase, which is 
consistent with the general expectation. For the two B — Kt]' decays, the large 
imaginary parts of TW"""* can also provide an effective enhancement to their branch- 
ing ratios. From this point we understand that although the factorizable emission 
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diagrams Figs. 1(a)- 1(b) provide the dominant contribution to the considered de- 
cays, but the LO contribution from the annihilation diagrams also provide essential 
contribution, therefore the NLO contribution Ai'^^Q from the corresponding anni- 
hilation diagrams as shown in Figs.5(e)-5(h) may be comparable with other NLO 
parts, and thus the analytical calculations for these diagrams should be done as 
soon as possible. 

Now we study the relative strength for the NLO contributions from different sources 
and collect all numerical results in Table IVIIII and IIXI 

In Table IVIIII we list the numerical values for individual decay amplitudes. The decay 
amplitude A4nlowc are obtained by evaluating the Figs. 1(a)- 1(h) with the use of the 
NLO Wilson coefficients Ci{fi), the NLO RG evolution matrix U{t,m,a) and as{fi) at 
two-loop level. The label Aivc denotes the changes of J^nlowc when only the NLO 
vertex corrections are included. The label Aiqi {J^mp ) shows the changes of J^nlowc 
when only the NLO contributions from the quark-loops ( the chromo-magnetic penguin) 
are included. The label Aipp shows the variation of M-nlowc when only the B ^ K 
and B — >■ r]^'^ transition form factors at NLO level are used. 

The label J^vc+qi+mp and M.nlo in Table IVnT] and J^nloi in Table HXl are defined 
as the following summation: 

M VC+ql+mp — Mvc + Mqi + Mmp, (89) 

■M.NLOI = -M-NLOWC + -M-YC+ql+mp-, (90) 

= Mnloi + Mff- (91) 
Where M.nloi is equivalent to the total decay amplitude M. as defined in Eq.(76) of 



Ref. |12|], and M.nlo is the decay amplitude when all currently known NLO contributions 
are taken into account. 

The ratio Rnlo in Table HXl is defined as 

Rnlo = \Mlo? ■■ \Mnloi? : l-M^^loP- (92) 



TABLE VIII. The numerical values (in unit of 10 ^) of the individual NLO contributions to the 
decay amplitudes, coming from different sources in the 77-77' and rj-rj'-G mixing scheme. 



Decays 


MS 


■M-NLOWC 


Mvc 




Mlmp 


MvC+ql+mp 


Mff 




1 

2 


-1.54 + i0.48 
-0.96 + i0.54 


-0.03 - iO.26 
-0.07 - iO.OS 


-0.32 - iOAA 
-0.25-i0.36 


0.35 - 70.27 
0.22 - 70.25 


-0.00-i0.97 
-0.10-i0.68 


0.05 - i0.04 
-0.001 - i0.02 


K%' 


1 

2 


4.83 - iSM 
5.62 - iAm 


0.63-il.34 
0.61 -il. 22 


1.29 + il.59 
1.40 + il.73 


-1.35 + i0.37 
-1.47 + i0.36 


0.57 + i0.62 
0.54 + i0.87 


0.50 - iO.25 
0.47 - iO.24 




1 

2 


-1.65 + iO.ll 
-1.13 + i0.06 


-0.05 - 70.35 
-0.10 - 7O.I8 


-0.32 - 70.44 
-0.25 - 70.36 


0.37 - 70.27 
0.26 - 70.26 


0.00-il.06 
-0.09-i0.80 


0.001 - 7O.I8 
-0.04-i0.13 


K+rj' 


1 

2 


4.66 - iAM 
5.39 - i4.57 


0.65-71.38 
0.62-71.29 


1.30 + il.59 
1.41 + il.73 


-1.38 + 70.36 
-1.50 + 70.35 


0.57 + i0.57 
0.53 + i0.79 


0.46 - 70.36 
0.43 - iO.35 



From the numerical results as shown in Tables IVIIIIIXI we find the following points: 
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TABLE IX. The numerical values of M.lo^ ■M.nloi-, J^nlo (in unit of 10 ^) and the ratio 
Rnlo for B — )■ Krj^'^ decays in the rj-Vj' and rj-rj'-G mixing schemes. 



Decay 


MS 


Mlo 


■M.NL01 


K/ltot 

i V L/Ly 


Rnlo 




1 

2 


-1.30 + i0.47 
-0.76 + i0.51 


-l.54-i0.49 
-1.05 - i0.14 


-1.49 - iO.53 
-1.05 - i0.15 


1.91 : 2.61 : 2.50 
0.84 : 1.12 : 1.13 


K^i 


1 
2 


3.40 - i3.29 
4.32 - i3.37 


5.41 -i3. 32 
6.16-^3.16 


5.91 - i3.57 
6.63 - i3.40 


22.4 : 40.3 : 46.3 
30.0 : 47.9 : 55.5 


K+7] 


1 
2 


-1.67- i0.06 
-1.05 - i0.03 


-1.65 - iO.96 
-1.22 - iO.74 


-1.65 - il.l3 
-1.26 - iO.87 


2.79 : 3.64 : 4.00 
1.10 : 2.04 : 2.35 


K+T]' 


1 

2 


3.50-i3.83 
4.31 - z4.08 


5.22 -i3.77 
5.91 -i3.78 


5.69 - i4.13 
6.34-i4.13 


26.9 : 41.5 : 49.4 
35.2 : 49.2 : 57.3 



(i) For all four B — )■ Kr]'^'^ decays, there are strong cancelation between M.vc^ -Mgi 
and A^mp- 

(ii) For two B — t- Kt] decays, the corresponding AipF are also much smaller in mag- 
nitude than other NLO parts M.vci -Mqi aiid Mmp, and also smaller in size than 
their summation Aivc+qi+mp- 

(iii) For two B — )■ Kr]' decays, the corresponding AipF are also much smaller than 
other NLO parts Mvc, M.qi and Aimp, but comparable in size with their summa- 
tion Aivc+qi+mp, and therefore all NLO contributions together provide the required 
enhancements to Br{B — )■ Krj') to account for the measured values. 

(iv) The only missing NLO parts in the pQCD approach are A^^^/q from Figs.5(a)-5(d) 
and Ai^^io from Figs.5(e)-5(h). They are most possibly small in size according to 
the studies in this paper and the general expectations based on the isospin cance- 
lation and power suppression. 

VI. SUMMARY 

In this paper, we made a systematic study about the four B Kt]^'^ decays in the 
pQCD factorization approach. We calculated the CP-averaged branching ratios and CP- 
violating asymmetries of the four B — i- Ki]^'^ decays in three different mixing schemes: 
the ordinary FKS t]--)]' mixing scheme, the r]-r}'-G and the rj-rj'-G-rjc mixing scheme. We 
considered the full LO contributions and all currently known NLO contributions to -B — )■ 
Krf^'^ decays in the pQCD approach. Besides those NLO contributions considered in 
Ref. [HI, we here also take the newly known NLO part of 5 — )■ (K, r^*^'^) transition form 
factors into account. 

From our numerical calculations and phenomenological analysis, we find the following 
points 

(1) In all the three considered mixing schemes, the NLO pQCD predictions for the 
branching ratios and CP-violating asymmetries agree with the data within one stan- 
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dard deviation, of course, partially due to the still large theoretical errors. But the 
NLO pQCD predictions in the tj-tj'-G mixing scheme provide a nearly perfect in- 
terpretation for the measured values. The NLO pQCD predictions in MS-2 are the 
following: 



Br{B^ K^7]') = {67.3tlli) x 10-^ (93) 



for branching ratios, and 



A'SpiB" KIt]') 



(-22.91^1?) X 10-^ 
(-5.5^}-^) X 10-^ 

(-16.0;if3) X lo-^ 
(66.7^1^) X 10-^ 
(3.3 ±0.3) X 10-^ 
(70.0 ± 0.3) X 10-^ 



for the CP- violating asymmetries, where the individual theoretical errors have been 
combined in quadrature. 

For B^ — )■ K^ri' and S+ — > K^ri' decays, the NLO contributions provide significant 
enhancements to their branching ratios. In the rj-rj'-G mixing scheme, for exam- 
ple, the NLO contribution provide a 89% (73%) enhancement to Br{B^ — )■ K^r]') 
{Br{B^ — )> K^rj') ) with respect to the LO prediction, such enhancements play the 
key role in our effort to resolve the i^?7*^'-*-puzzle and to understand the patten of 
the Br{B Krj^'^). 

For B^ — > K^T] and B^ — )> K^rj decays, although the inclusion of the NLO contri- 
butions only leads to a relatively small changes to their branching ratios, but the 
resulted variations are in the right direction and helpful for us to improve the consis- 
tency between the pQCD predictions and the measured values. In the rj-rj'-G mixing 
scheme, for instance, the pQCD predictions are Br{B^ — )■ K^r]) = 0.90 x 10^^ and 
Br{B^ K^rf) = 1.98 x 10^^ at the leading order, and changed into Br{B^ 
K^i) = 1.13 X 10-^ and Br{B+ K+rj) = 2.36 x 10"^ when the NLO contributions 
are taken into account, while the corresponding measured values are 1.23^q|23 xlO-^ 
and 2.36lg|? x 10"^ respectively. 

By comparing the pQCD predictions as given in "NLOl" and "NLO" columns in 
Tables II-IV, one can see directly the effects of NLO form factors: the NLO part 
Ai FF oi the B ^ K and B — )■ 77^'^ form factors can produce about 20% enhancement 
to the branching ratios Br{B — )■ Kt]'), which plays an important role in closing the 
gap between the pQCD predictions and the relevant data. 
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(5) In the r]-r]'-G-r]c mixing scheme, the decay chain B — )■ Kr]c -n- Kr]^'^ can provide 
an effective enhancement to the branching ratios at the leading-order. But when 
the large NLO contributions are taken into account, the effects of the ric component 
become unimportant. 

(6) For — ^ K^T] decays, the LO pQCD predictions for in all three mixing 
schemes have an opposite sign with the measured value. The inclusion of the NLO 
contributions changed the sign of the pQCD predictions for A^p, while the NLO 
pQCD predictions for A'cp{B^toK^r]) in the cases of MS-1 and MS-2 are now 
becoming consistent with the data within one standard deviation, but the NLO 
pQCD prediction for ^^p(i?^ — > K^rf) in the case of MS-3 is still much smaller 
than the measured value in magnitude. 

(7) For A'cp{B^toK^r}'), the NLO pQCD predictions agree with the data within one 
standard deviation, while the consistency between the pQCD predictions and the 
data are improved by the inclusion of the NLO contributions. 

(8) For the direct and mixing-induced CP-violating asymmetries Ac}>{B'^ — >■ Kgi]^'^) 
and Acp{B° — > Kgrj^'^), the pQCD predictions have weak dependence on the 
NLO contributions and the choice of different mixing schemes. For Aqp"^^^{B'^ — )■ 
KgT]'), for example, the NLO pQCD predictions are Acp{B° ^sV') ~ 3% and 
Aqp{B^ — )■ K^Tj') ~ 70%, which are well consistent with the measured values of 
(1 ± 9)% and (64 ± 11)% respectively. 

(9) The factorizablc emission diagrams Figs.l(a)-l(b) provide the dominant contribu- 
tion to the considered decays. The LO contribution Ai'^^'^ from the spectator dia- 
grams Figs.l(c)-l(d) is already less than 3% of the total contribution, the next-to- 
leading order contributions M.'^^^q from the diagrams Figs.5(a)-5(d) are the higher 
order contribution and therefore should be smaller than its LO counterpart. Con- 
sequently, it is reasonable for us to neglect M.'^^lo from the spectator diagrams 
Figs.5(a)-5(d). 

(10) The real part of is always neghgibly small, but its imaginary part is relatively 

large and leads to a large strong phase, which can also produce an effective enhance- 
ment to the branching ratios of the considered decays. Although |A^^£ol most 
possibly much smaller than its LO counterpart |A^""™|, but the still missing NLO 
contribution A4'}^lo i^om Figs.5(c)-5(li) may be comparable in size with A4ff, and 
should be calculated as soon as possible. 
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Appendix A: Distribution Amplitudes 



The expressions of the relevant distribution amphtudes (DA's) of K meson are the 
following (33[: 



^-^6x(l -x)[l + af + af (t) + af (t) 

fK 



(Al) 



9 



1 + 6(5r/3 - Imu^s - ^Pk - ^Pi-af )(1 - 10a; + lOx^) 



(A3) 



with the mass ratio px = friK/'^oK- The Gegenbauer moments are of the form [33 

-0.015. 



af = 0.2, af = 0.25, af 



The values of other parameters are rj^ = 0.015 and u 
polynomials C^(t) are given as: 



(A4) 

3.0. At last the Gegenbauer 



C^/\t) = -{3t' - 1), C^^\t) = -(3 - 30t^ + 35t 



1/2, 



Cf (t) = 3t, C','\t) = - 1), 



.3/2, 



15 



:i - i4r + 2ir 



(A5) 



with t = 2a; — 1. 



The distribution amplitudes (p^f'^ are given as [33 



2v^ 

+al^Cl'\2x-l) 



6x(l - x) 1 + a^'Ci/^(2x - 1) + al'C^^'\2x - 1) 



2v/2iVc 



1 + (30r/3 - ^pyC2^/'(2a; - 1) 
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-3 <! r?3C^3 + ^Pi(l + 6a^') |> Cl/'(2a; - 1) 



;i - 2x) 



■ (1 - 10x + lOx^)] 



1 7 3 

1 + 6 ( 5r/3 - 2 ^3^3 - — pj, - -p^al' 



(A6) 



(A7) 



(A^ 



Vq 



0, a, 



0.44 ± 0.22, a 



Vq 
, 014 



= 0.25, and the 



where p^^ = 2mg/mgg,al 
Gegenbauer polynomials C^{t) have been given in Eq. flA5p . As to the wave function and 
the corresponding DA's of the ss components, we also use the same form as qq but with 
some parameters changed: p^^ = 2ms/mss, aT = al'' for i = 1,2,4. 
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Appendix B: Related Hard Functions 

In this appendix, we present the expUcit expressions for those hard scales and hard 
functions. 

The hard scales appeared in the decay amplitudes are chosen as 

ta = max{y/x^MB, l/h, 1/h}, 
t'^ = max{y/xlMB, 1/bi, 1/63}, 

tb = max{y/xiX3MB, a/|1 - Xi - X2\x3Mb, l/h, I/62}, 
t[ = max{y/x^MB, >/\x^~^^x^\x^Mb, l/h, I/62}, 
tc = max{Vl - XsMb, I/&2, 1/&3}, 
t'^ = max{^/x^MB, I/&2, ^/h}, 



td = max{A/.T2(l - Xs)Mb, ^/l - {1 - Xi - X2)x3Mb, l/^i, I/&2}, 
t'^ = max{^/x^{T^^X3)MB, ^/\x^~^^x^\(T^^XsjMB, 1/61, 1/62}, 



te = max{^,T3(l - r2 )Mb, 1/h, 

t'^ = max{^xi(l-r2jMs, 1/h: V&s}, 

tf = max{^xiX3{l - t^JMb, ^\{-l + Xi +a;2)[a;3 + {1 - X2 - X3)r^J +r^jMB, 

l/h,l/b2}, 

t'f = max{y^xiX3(l - r^jM^, ^J\{xi- X2)[x3 + {x2 - X3)r^J + v^JMb, 

l/&i,l/fc2}. (Bl) 

The hard functions hiS appeared in the decay amplitudes are defined by 

he{xi,X3,h:h) = [0{h - h) loiVx^Meh) Koi^/^Meh) 

+e(h - h)Io{V^MBh)Ko{^MBh)] Ko(^/^MBh)St(x3), 
h„{xi,X2,X3,h,h) = [d{h - h)Ko{^xiX3MBh)h{\/xiX3MBh) 

+e(h - h)Ko(^/^MBh)Io(V^MBh)] 

^ ■ fH^^\^/{x2 - Xi)x3MBh), Xi-X2<0 ^g2) 

Ko{-\/ (Xi - X2)x3MBh), Xi-X2>0 



ATT, 



ha{x2,X3,h,h)^CffSt{x3) 9 {h - h) H^o\V^^M sh) UV^^M sh) 



+eih - h)H^'\V^MBh)MV^MBh) H'^''iy/^MBh), 



r(i)/ 



hna{Xl,X2,X3,h,h) = y 0{h - h)H^^\\/x2{l - X3)MBh)MV^2{l - X3)MBh) 

+9{h - h)H^^\VM^ - X3)MBh)MVM^ - X3)MBh) 

xi^o( ~{l-xi- X2)x3MBbi), (B3) 
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h'^^{xi,X2,X3,bi,b2) = — 0{bi - b2)H^o^\^/x2{l - X3)Mb&i) 70(^3^2(1 - 3^3)^562 

r(i) 



+^(62 - b^)Hi,'\^/x2il - x,)MBb2)MV^2il - X:,)MBbi] 

f g^' V(x2 - Xi)(l - X3)MBbi), Xi - X2 < 
-^o(a/ {Xi - X2){1 - X3)Mb6i), Xi - X2 > 



(B4) 



where = h{z) + iYo{z). 

K{xuXs,biM) = [e{bi-b3)Io{y^MBb3)Ko{V^MBbi) 

+e{bs - b,)h{V^MBbi)Ko{^MBbs)] Ko{,/pMBbi)St{xs), 



(B5) 



where a = Xs{l - r^J, /3 = XiXs^l - r^J. 

/i;(a;i,X2,X3,6i,62) = [^(62 - &i)^o(v^Mb62)/o(v^Mb6i) 

+^(61 - b2)Ko{V^MBbi)Io{V^MBb2) 



X 



Ko{,/\fi\MBb2), 0^ > 

where a' = xia;3(l — r^J, = (xi — X2)[x2r'^^ + xifl — r^J] + r^^. 
The function S't(a;) has been parameterized as (34| 

2i+2T(3/2 + c) 



(B6) 



(B7) 



v^r(i + c 

with c = 0.3. 

The evolution factors Ee \ and i?a\ appeared in the decay amphtudes are given by 

Ee{t) = «,(t)exp[-^B(t) - S^{t)l 

E'^it) = a,{t)exp[-SB{t) - S2{t) - S3{t)]\b,=bs, 

Eait) = a,{t)exp[-S2{t) - S^it)], 

E'^it) = as{t)eM-SB{t) - S2{t) - S3{t)]\b,=H. (B8 
where the Sudakov exponents are defined as 



SB{t) = s\xi 

S2{t) = S I X2 



Mb 

Mb 
V2' 



— 7g(as(/^)), 



<J Jl/bi 
+ S ( (1 - X2)^,&2 



V2' 



+ 2 



— 7g(«s(/i)), 

1/62 



(B9) 



with the quark anomalous dimension 7^ = —as/ir. Replacing the kinematic variables of 
M2 to M3 in 5*2, we can get the expression for 5*3. The explicit form for the functions 



s{Q,b) can be found in Ref. [21 
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